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Abstract. Generically, in models of driven interacting particles the average 
speed of the particles decreases monotonically with increasing density. We propose 
a counter-example, motivated by the motion of ants in a trail, where the average 
speed of the particles varies non-monotonically with their density because of the 
coupling of their dynamics with another dynamical variable. These results, in 
principle, can be tested experimentally. 

Particle-hopping models, formulated usually in terms of cellular automata (CA) 
, have been used to study the spatio-temporal organization in systems of interacting 
particles driven far from equilibrium [2 03 El 03 which include, for example, vehicular 
traffic Ej • In general, the inter-particle interactions tend to hinder their motions 
so that the average speed decreases monotonically with the increasing density of the 
particles. In this letter we report a counter-example, motivated by the motion of 
ants in a trail [5], where the average speed of the particles varies non-monotonically 
with their density because of the coupling of their dynamics with another dynamical 
variable. 

The ants communicate with each other through a process called chemotaxis by 
dropping a chemical (generically called pheromone) on the substrate as they crawl 
forward E3- Although we cannot smell it, the trail pheromone sticks to the 
substrate long enough for the other following sniffing ants to pick up its smell and 
follow the trail. In this letter we develope a CA model which may be interpreted as 
a model of unidirectional flow in an ant-trail. Rather than addressing the question of 
the emergence of the ant-trail, we focus attention here on the traffic of ants on a trail 
which has already been formed. 

Each site of our one-dimensional ant-trail model represents a cell that can 
accomodate at most one ant at a time (see fignj. The lattice sites are labelled 
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Figure 1. Schematic representation of typical configurations; it also illustrates 
the update procedure. Top: Configuration at time t, i.e. before stage / of the 
update. The non-vanishing hopping probabilities of the ants are also shown 
explicitly. Middle: Configuration after one possible realisation of stage I. Two 
ants have moved compared to the top part of the figure. Also indicated are 
the pheromones that may evaporate in stage // of the update scheme. Bottom: 
Configuration after one possible realization of stage II. Two pheromones have 
evaporated and one pheromone has been created due to the motion of an ant. 



by the index i (i = 1,2, ...,L); L being the length of the lattice. We associate 
two binary variables Si and Uj with each site i where Si takes the value or 1 
depending on whether the cell is empty or occupied by an ant. Similarly, ai = 1 
if the cell i contains pheromone; otherwise, o~i = 0. Thus, we have two subsets of 
dynamical variables in this model, namely, {S(t)} = (S\(t), S2(t), Si(t), SL(i)) 
and {o-(t)} = (<ri(t), <72(t), cn(t), 0z,(f)). The instantaneous state (i.e., the 
configuration) of the system at any time is specified completely by the set ({5*}, {<x}). 

We assume that the ant does not move backward; its forward-hopping probabil- 
ity, however, is higher if it smells pheromone ahead of it. The state of the system is 
updated at each time step in two stages. At the end of stage I we obtain the subset 
{S(t + 1)} at the time step t + 1 using the full information ({S(t)}, {cr(t)}) at time t. 
At the end of the stage II we obtain the subset {a(t + 1)} at the time step t + 1 using 
the subsets {S(t + 1)} and {a(i}}. 

Stage I: The subset {5*} (i.e., the positions of the ants) is updated in parallel according 
to the following rules: 

If Si(t) = 1, i.e., the cell i is occupied by an ant at the time step t, then the ant hops 
forward to the next cell i + 1 with 

r Q if S i+ i(t)=0butCT i+ i(t) = l, 
probability = < q if Sj+i (t) = and a-j+i (t) = 0, (1) 
{ if S i+1 (t) = l. 
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where, to be consistent with real ant-trails, we assume q < Q. 

Stage II: The subset {a} (i.e., the presence or absence of pheromones) is updated in 
parallel according to the following rules: 

If <7j(f) = 1, i.e., the cell i contains pheromone at the time step t, then it contains 
pheromone also in the next time step, i.e., aiit +!) = !, with 



where / is the pheromone evaporation probability per unit time. 

On the other hand, if o~i(t) — 0, i.e., the cell i does not contain pheromone at the time 

step t, then 



In certain limits our model reduces to the Nagel-Schreckenberg (NS) model 11 
which is the minimal particle-hopping model for vehicular traffic on freeways. The 
most important quantity of interest in the context of flow properties of the traffic 
models is the fundamental diagram, i.e., the flux-versus-density relation, where flux is 
the product of the density and the average speed. For a hopping probability q^s at a 
given density c the exact flux F(c) in the NS model is given by |51ll2j 



which reduces to -Fjvs(c) = min(c, 1 — c) in the deterministic limit q^s = !• 

Note that in the two special cases / — and / = 1 the ant-trail model becomes 
identical to the NS model with qpjs — Q and q^s = Q, respectively. Extensions 
of the NS model have been used not only to capture different aspects of vehicular 
traffic [Hj but also to simulate pedestrian dynamics In a closely related CA 

model for pedestrian dynamics Jl] the floor fields, albeit virtual, are analogs of the 
pheromone fields {a} in the ant-trail model. However, in the pedestrian model there 
is no exclusion principle for the floor field. 

The ant-trail model we propose here is also closely related to the bus route model 
(BRM) |15l ni)| . The variables S and a in the ant-trail model are the analogs of the 
variables representing the presence (or absence) of bus and passengers, respectively, 
in the BRM. Because of the periodic boundary conditions, the number of ants and 
buses are conserved while the pheromone and passengers are not conserved. However, 
unlike the BRM, the pheromones are not dropped independently from outside, but 
by the ants themselves. Another crucial difference between these two models is that 
in the bus-route model Q < q (as the buses must slow down to pick up the waiting 
passengers) whereas in our ant-trail model Q > q (because an ant is more likely to 
move forward if it smells pheromone ahead of it). 

In figEl we show the fundamental diagrams obtained by extensive computer 
simulations of the ant-trail model for several values of /. The most unusual features 
of the fundamental diagrams shown in figElare that, over an intermediate range of 

% By the term 'NS model' in this letter we shall always mean the NS model with maximum allowed 
speed unity, so that each particle can move forward, by one lattice spacing, with probability qjvs if 
the lattice site immediately in front is empty. 




<Ti(t + 1) = 1 iff Si(t + 1) = 1 at the end of stage I. 



(3) 



Fns{c) = - [l - V 1 - 4 Ins c(l - c) 



(4) 
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Figure 2. The flux of ants plotted against their densities for the parameters 
(a) Q = l,q = 0.25 and (b) Q = 0.75, q = 0.25. The discrete data points 
corresponding to / = 0.0001(v), 0.0005(0), 0.001(o), 0.005(«), 0.01(A), 0.05(D), 
0.10(x), 0.25(+), 0.50(*) have been obtained from computer simulations; the 
dotted lines connecting these data points merely serve as the guide to the eye. 
The two continuous solid curves at the top and bottom correspond to the flux 
in the NS model for q^s = 1-0 and q^s = 0.25, respectively, in (a) and for 
qpfg = 0.75 and qpfs = 0.25, respectively, in (b). 

values of / (for example, / = 0.0005,0.001,0.005,0.01 in figgjl the flux in the low- 
density limit c — > is very close to that for the NS model with q^s = q whereas in 
the high-density limit c — > 1 the flux for the same / is almost identical to that for the 
NS model with qns = Q- These unusual features of the fundamental diagrams arise 
from the non-monotonic variation of the average velocity with the density of the ants 
(see fig EH). 

The presence of the pheromone essentially introduces an effective hopping 
probability g e g(c), which depends on the ant density c. The particle-hole symmetry 
(and hence the symmetry of the fundamental diagram about c = 1/2) observed in 
the special limits / = and / = 1, are broken by the c-dependent effective hopping 
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Figure 3. The average velocity of ants plotted against their densities for the 
parameters (a) Q = 1, q = 0.25 and (b) Q = 0.75, q = 0.25. Same symbols in the 
figs|5] and |21 correspond to the same values of the parameter /. 



probability for all < / < 1 leading to a peak at c > 1/2. Furthermore, the analysis 
of correlation functions reveals some interesting clustering properties which will be 
studied in detail in a future publication |17| . 

The qualitative features of the c-dependence of g e ff can be reproduced by an 
analytical argument based on a mean-field approximation (MFA) 17 . In this MFA, 
let us assume that all the ants move with the mean velocity (V) which depends on 
the density c of the ants as well as on /; although, to begin with, the nature of these 
dependences are not known we'll obtain these self-consistently. 

Let us consider a pair of ants having a gap of n sites in between. The 
probability that the site immediately in front of the following ant contains phcromonc 
is (1 — f) n /( v > since -jSp- is the average time since the pheromone has been dropped. 
Therefore, in the MFA the effective hopping probability is given by 

fcff = Q(l - I) n/(V) + Q{1 - (1 - f) n/(V) }. (5) 
We replace n by the corresponding exact global mean separation (n) = - — 1 between 
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successive ants. Moreover, since (V) is identical to q D g, we get the equation 




(6) 



which is to be solved self-consistently for q e g as a function of c for a given /. Note 
that the equation © implies that, for given /, lim c ^o q e g = q; this reflects the fact 
that, in the low-density regime, the pheromone dropped by an ant gets enough time to 
completely evaporate before the following ant comes close enough to smell it. Equation 
JfjJ) also implies lim c ^i q c g — Q; this captures the sufficiently high density situations 
where the ants are too close to miss the smell of the pheromone dropped by the leading 
ant unless the pheromone evaporation probability / is very high. Similarly, from © 
we get, for given c, lim/_,i q e g = q and lim/_,o QeS = Q which are also consistent with 
intuitive expectations. 

In view of the fact J0| that the lifetime of pheromones can be as long as thirty 
to sixty minutes, the interesting regime of / 1), where the average velocity varies 
non-monotonically with the ant density, seems to be experimentally accessible. We 
hope that the non-trivial predictions of this minimal model of ant-trail will stimulate 
experimental measurement of the ant flux as a function of the ant density for different 
rates of pheromone evaporation by using different varieties of ants |S] . 
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